How sensitive is a neutrino factory to the angle #13? 
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We consider the impact of non-standard interactions of neutrinos (NSI) on the determination 
of neutrino mixing parameters at a neutrino factory using v e — > "golden channels" for the 
measurement of #13. We show how a small residual NSI leads to a drastic loss in sensitivity in #13, 
of up to two orders of magnitude. This can be somewhat overcome if two baselines are combined. 



The long-standing solar and atmospheric neu- 
trino anomalies Jl| now give compelling evidence 
for new physics in the lepton sector. The data are 
well described in the simplest 3-neutrino oscillation 
framework Q defined by two mass-squared differ- 
ences Am 2 ol and Am 2 tm , three angles #12, 623, 6*13 
and one relevant CP-violating phase S. The pa- 
rameters Am 2 tm and #23 are determined by at- 
mospheric experiments, while Am 2 ol and #12 are 
determined by solar experiments. For $13 there is 
only an upper bound sin 2 2^13 < 0.1 at 90% CL 
derived from reactor experiments H . The phase 5 
is completely unknown. 

A new generation of long-baseline neutrino oscil- 
lation experiments using a neutrino beam from the 
decay of muons in a storage ring is discussed 
Such so-called neutrino factories should give 
a precise determination of #13, possibly down to 
fewxlO -4 , exploring also the important issue of 
leptonic CP violation ||. 

In a large class of models neutrino masses are 
accompanied by non-standard interactions of neu- 
trinos (NSI). These can be non- universal (NU) or 
flavour-changing (FCI), and may arise also in the 
absence of neutrino mass |^|. Such non-standard 
interactions of neutrinos affect their propagation 
in matter, with implications for solar (7[ ||, ||, pxj ], 
atmospheric JO, and astrophysical neutrino 
sources Q. The possibility to use a neutrino 
factory for probing NSI with the earth matter, 
which neutrinos have to cross necessarily in a long- 
baseline experiment, has been previously examined 
in the appearance channels v e — ► v T [H and 

In this letter we will consider the impact of NSI 
on the determination of neutrino mixing param- 
eters assuming the standard setup of a neutrino 
factory. Therefore, we will focus on the v e — ► ( z/^ 
channels, i. e. the "golden channels" for the mea- 
surement of #13 . In contrast to refs. jl4|, we 



need only to assume a muon detector, which seems 
more realistic than a tau detector. We treat NSI 
and neutrino oscillations in a common framework 
to investigate correlations between the correspond- 
ing parameters. We show that the presence of a 
small FCI in the e — r channel, safely within cur- 
rent limits, leads to a drastic loss in sensitivity in 
013, up to two orders of magnitude. This follows 
from a strong correlation between these parame- 
ters. We show how this effect can be diminished 
to some extent if two baselines are combined. 

In the presence of NSI the evolution of neutri- 
nos passing through matter is modified by a non- 
standard potential due to the coherent forward 
scattering amplitude of NSI-processes v a + / — > 
v fi + /1 where / is a fermion in the medium. For 
simplicity we will consider only NSI with down- 
quarks: / = d. Then the potential induced by NSI 
can be written as Unsi = V%Gf Nd t a p, where Nd 
is the number density of down-quarks along the 
neutrino path. The small numbers e a p describe 
the NSI in units of the Fermi-constant Gf- The 
off-diagonal elements e a p with a ^ (3 correspond 
to FCI, whereas the differences in the diagonal el- 
ements e aa lead to NU interactions. For simplic- 
ity we take all t a p real. Present data on various 
lepton flavour violating processes can be used to 
set upper bounds on FCI |l^, Il6| [TtJ leading to 
M <7xl0- 5 , | C|iT | ^xlO-rfcer! < 7xHT 2 . 
Much weaker bounds of order 0.1 apply for the NU 
coefficients. 

In this letter we show how claimed sensitivities 
on 613 get destroyed if a small NSI component ex- 
ists. For this it suffices to consider a simplified sce- 
nario in which we neglect the solar mass-squared 
difference. If we parameterize the lepton mixing 
matrix U23U13U12, this implies also that the angle 
0\2 and the phase S disappear. Further, we restrict 
ourselves to the elements e eT and e TT (the case e^ T 
and e Tr is considered in Ref. |l5f| ), where present 
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constraints are rather weak. Then the evolution of 
flavour neutrinos is governed by the Hamiltonian 

H u = U 23 U 13 di&g (0,0, A) U} 3 U1 3 
1 re eT \ 
-V I . (1) 
rt eT re TT J 

Here we have defined A := Am^ to /2£„ and 
V = ^/2GfN c is the matter potential due to the 
Standard Model charged current interaction ||], 
where N e is the electron number density and r := 
Vnsi/V = Nd/N e with r « 3 in earth matter. In 
the Hamiltonian ([!]) a sign change of A is equiva- 
lent to a sign change of V, which interchanges the 
evolution of neutrinos and anti-neutrinos. In the 
following we show that even in this simple scenario 
the strong correlation between S13 and e eT will lead 
to serious complications for the measurement of # 13 
at a neutrino factory. 

To get an intuitive understanding of the numer- 
ical results it is useful to derive first an analytic 
expression for the appearance probability. To this 
aim we assume a constant matter potential V and, 
since S13 must be small, consider all terms contain- 
ing S13, e eT and e TT as a small perturbation of the 
Hamiltonian and calculate eigenvalues and eigen- 
vectors of eq. [j] up to first order in these small 
quantities. The resulting appearance probability 
is second order and has the form 

P„ e ^ w As 2 13 + Bs 13 e eT + Ce 2 eT . (2) 

where the coefficients A, B and C are of the same 
order and depend on the neutrino energy, on the 
baseline L and on the sign of V (neutrinos or anti- 
neutrinos) in a nontrivial way. For example, 
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with similar expressions for B and C . These an- 
alytic expressions for the appearance probabili- 
ties are in agreement with numerical calculations 
within a few % in the relevant parameter range. 
We note that the NU coefficient e TT does not show 
up in leading order in the appearance probability 
[ p"4| . Therefore, we will neglect it in our numerical 
studies and concentrate only on FCI. We consider 
both signs for e eT , this sign however becomes irrel- 
evant as S13 approaches zero. 

Eq. ^| implies that the lines of constant event 
rates in the Si3-e eT plane are ellipses, whose coef- 
ficients are obtained by fixing a baseline and fold- 
ing A, B and C with flux and detection cross sec- 
tion. In fig. [l] we draw these ellipses for neutrinos 




Figure 1: Lines of constant event rates in si3-e eT 
(black for neutrinos and grey for anti-neutrinos) at 
L = 3 000 km. 



and anti-neutrinos through different start points 
marked with a star. In general the two ellipses 
will cross not only at the start point. All points 
where the ellipses cross will lead to identical total 
event rates and are therefore a source of possible 
confusion. Depending on the choice of initial val- 
ues four typical cases can be distinguished. Case 
(a) has two solutions witch are very close and away 
from both axes. In this case both NSI and oscil- 
lations can be established. Case (b) has a second 
solution on the e eT = axes. Therefore one can 
only establish the existence of oscillations. Case 
(c) has additional solutions on the S13 = axes, 
which means that only NSI can be established but 
not oscillation. Case (d) has solutions which are 
close to both axes. If the count rate is not high 
enough the allowed region will touch both axes, 
leading to complete confusion between NSI and os- 
cillation. In the following we explore this through a 
detailed numerical example and show that all these 
four cases can actually arise for realistic oscillation 
and NSI parameter choices. 

In our numerical calculations we assume a neu- 
trino factory with an energy of 50 GeV for the 
stored muons and 2 x 10 20 useful muon decays 
of each polarity. We consider a magnetized iron 
detector with perfect charge identification of the 
muons and a mass of 10 kt. The muon detection 
threshold is set to 4 GeV, the energy resolution of 
the detector is approximated by a Gaussian reso- 
lution function with AE V /E U — 10% and we use 
20 bins in muon energy. The appearance probabil- 
ity is obtained by numerically solving the neutrino 
evolution equation with the Hamiltonian (El), using 
a realistic earth matter density profile. Then this 
probability is folded with neutrino flux, cross sec- 
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tion and energy resolution function to obtain the 
event rates in the detector, following |f[gj| . 

We consider only the v e — > appearance chan- 
nel since it contains most of the information. Our 
"observables" are the event rates of this channel 
n% u i n h) f° r (anti-)neutrinos in each energy bin i. 
We fix the atmospheric oscillation parameters at 
Am 2 tm = 3 x lCT 3 eV 2 and sin 2 26> 23 = 1. Thus, at 
a given baseline the event rates depend on S13 := 
sin 2 2^i3 and e eT . In order to evaluate the poten- 
tial of a neutrino factory to probe these parameters 
we proceed as follows. For any given (5° 3 , e° T ) 
choice we calculate the event rates n l x (Si 3 , e° eT ) 
with x = v, v and construct a X (<Si3i e ei-; S® 3 , £g T ) 
appropriate for a Poisson distribution. Using this 
we derive allowed regions in the Si 3 — e eT plane as 
usual. For all S , 13 and e eT values inside these con- 
fidence regions the data can be fitted at the given 
CL (in our case 99%) if the true values are Sf 3 and 
£g T . Thus, these regions indicate the capability of 
a neutrino factory to probe Si 3 and e eT . 

In fig. U we show such fits for four characteristic 
test points (S° 3 ,eg r ) at a baseline of 3 000km in 
correspondence with the cases considered in fig. |l|. 
The two grey lines passing through each of these 
test points are the lines of constant total rates for 
neutrinos and anti-neutrinos. From the fact that 
the confidence regions follow closely the lines of 
constant rates we learn that most information is 
contained in the total rates; the spectral informa- 
tion is not very important. We have verified that 
our results are, indeed, rather insensitive to varia- 
tions of the number of energy bins and of the en- 
ergy resolution. Moreover we find that one must 
take into account simultaneously both neutrino 
and anti-neutrino rates (this requires to run the 
neutrino factory in both polarities). This is easy 
to understand since the allowed regions extend as 
long as the lines of constant neutrino and anti- 
neutrino rates remain close to each other, i. e. both, 
rates are similar to the ones at the test point. In 
regions where the lines are far from each other the 
neutrino or the anti-neutrino rates are very differ- 
ent from the one in the test point, and hence it is 
not possible to fit the data generated by the pa- 
rameters (S° 3 , e° T ) in those regions. 

Let us consider in more detail the fits for the 
four different test points in fig. |[ For case (a) with 
large S° 3 and e° T values one would need both os- 
cillations and FCI. In case (b) (S? 3 = lO" 2 , e° r = 
10~ 2 ) the allowed region extends to the e eT = 
axis. This means that the data in a neutrino fac- 
tory generated by these parameters require oscilla- 
tions. In contrast, for case (c) (S° 3 = I0~ 3 , e a eT — 
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cated by the points, at a baseline of 3 000 km. Lines of 
constant event rates are displayed in grey. 



3.2 x I0 -2 ) the allowed region extends only to the 
Si 3 = axis, indicating that the data in a neu- 
trino factory generated by these parameters re- 
quire the presence of FCI. Finally, in case (d) 
(5° 3 = 10~ 3 , e" T = I0~ 2 ) there is complete con- 
fusion between oscillations and FCI, because the 
allowed region extends to both axes e eT = and 
Si 3 = 0. This means that the data in a neutrino 
factory generated by these parameters can be fit- 
ted either with pure oscillations (e eT = 0) or pure 
FCI (S13 = 0) with both signs of e eT . We can 
define a sensitivity to S13 by requiring that the al- 
lowed region for a given point in the S13 — e eT plane 
does not touch the Si 3 = axis. By scanning the 
S13 — e eT plane we obtain the following sensitivity 
limits for three baselines: 



baseline 


700 


3000 


7000 


w/o NSI 


2 x I0~ 4 


3 x I0" 4 


5 x I0" 4 


w NSI 


I x 10- 1 


7 x I0~ 2 


5 x I0~ 3 



Comparing with the sensitivity limits without 
NSI W we can see that the sensitivity loss is dra- 
matic, especially for shorter baselines. 

The geometry of the muon storage rings cur- 
rently discussed offers the striking possibility to 
illuminate two detectors at different baselines with 
neutrinos from one neutrino factory. In fig. 3 
we show the 99% CL region from a fit to data 
generated from the test point SJ3 = 10~ 3 and 
e° r = 10~ 3 at 700 km, at 7000 km and from a com- 
bination of both. In this particular case, for each 
of the two baselines taken alone there is complete 
confusion between S13 and e eT . However, combin- 
ing the two baselines we obtain a much smaller 
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work |l4| stimulated our interest. 
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Figure 3: sin 2 2#i3-e eT sensitivities when combining 
two baselines at 700 km, at 7000 km, see text. 



allowed region indicated by the solid line, which is 
roughly given by the intersection of the individual 
regions. Here the data can be fitted with e er = 0, 
with improved sensitivity for 6*13. For the combi- 
nation of two baselines the analysis is performed as 
before, but using the sum of the % 2 -functions for 
the individual baselines. We have evaluated neu- 
trino factory sensitivities (defined as above) for all 
possible combinations of the three baselines 700 
km, 3000 km and 7000 km. 



baseline 



w/o NSI 
wNSI 



700&3000 700&7000 3000&7000 



2 x 10~ 4 3 x 10~ 4 3 x 10~ 4 
2 x 10~ 3 2 x 10~ 3 2 x 10~ 3 



The sensitivities for S13 are still considerably worse 
(up to one order of magnitude) than those for the 
pure oscillation case (no NSI). Notice that no sub- 
stantial improvement is expected by changing the 
energy of the stored muons since neutrinos carry a 
relatively wide spectrum already for our nominal 
50 GeV choice. 

We have shown within a simplified scenario that 
the presence of even a small NSI component, pre- 
dicted in most models of neutrino mass, can deteri- 
orate the sensitivity of a neutrino factory for #13 by 
orders of magnitude. This problem can partly be 
resolved by combining two baselines. Additional 
effects of NSI in neutrino source and detection Q , 
and CP violating phases in NSI coefficients, as well 
as in the mixing matrix will lead to even more se- 
rious difficulties p0| . 
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